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Abstract 

In this article we formulate and discuss one particle quantum scattering theory on an ar- 
bitrary finite graph with n open ends and where we define the Hamiltonian to be (minus) the 
Laplace operator with general boundary conditions at the vertices. This results in a scatter- 
ing theory with n channels. The corresponding on-shell S-matrix formed by the reflection 
and transmission amplitudes for incoming plane waves of energy £' > is explicitly given 
in terms of the boundary conditions and the lengths of the internal lines. It is shown to 
be unitary, which may be viewed as the quantum version of Kirchhoff's law. We exhibit 
covariance and symmetry properties. It is symmetric if the boundary conditions are real. 
Also there is a duality transformation on the set of boundary conditions and the lengths of 
the internal lines such that the low energy behaviour of one theory gives the high energy 
behaviour of the transformed theory. Finally we provide a composition rule by which the 
on-shell S-matrix of a graph is factorizable in terms of the S-matrices of its subgraphs. 
All proofs only use known facts from the theory of self-adjoint extensions, standard linear 
algebra, complex function theory and elementary arguments from the theory of Hermitian 
symplectic forms. 
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1 Introduction 



At present mesocopic quasi-one-dimensional structures like quantum ||T], ^, atomic [^] and 
molecular |Q] wires have become the subject of intensive experimental and theoretical studies. 
This kind of electronics is still far from being commercially useful. However, the enormous 
progress that has been made in the past years suggests that it will not be too long before the first 
molecule-sized electronic components become a reality (see e.g. §,01). 

According to already traditional physical terminology a quantum wire is a graph-like struc- 
ture on a surface of a semiconductor, which confines an electron to potential grooves of width of 
about a few nanometers. An accurate theory for these nanostructures must include confinement, 
coupling between closely spaced wires, rough boundaries, impurities, etc. The simplest model 
describing the conduction in quantum wires is a Hamiltonian on a planar graph. A similar model 
can be applied to molecular wire - a quasi-one-dimensional molecule that can transport charge 
carriers (electrons or holes) between its ends [||]. Atomic wires, i.e. lines of metal atoms on 
the surface of a semiconductor provide another example of such quasi-one-dimensional struc- 
tures. Also Hamiltonians on planar graphs arise naturally in the modelling of high-temperature 
granular superconductors [Q, |l^, [ll]]. 

Although such models were proposed long ago (see e.g. [12, 13]), probably it was Pavlov 
and Gerasimenko [14, 15 1 who initiated a rigorous mathematical analysis of such models, which 
later acquired the name of quantum wires. A more general approach to the problem of the cor- 
responding mathematical structure was formulated in [|l^] several decades earlier. Here we do 
not intend to give a complete overview of the whole subject. We only mention some related 
studies. 



In [|17|] and [ |18| ] networks with leads were used to study adiabatic transport and Chern 
numbers. Two particle scattering theory on graphs was studied in [Jl9|]. Quantum waveguides 
[2C, 21, 2^, 24, 25, 2^, where the influence of confining potentials walls is modeled by the 
Dirichlet boundary condition, give a more realistic description of quasi-one dimensional con- 
ductors. The wave function is allowed to have several mutually interacting transversal modes. 
In real quantum wires the number of these modes can be rather large (up to 10^ - 10^). For 
another more realistic model of a two dimensional quantum wire see [27]. 

In this article we consider idealized quantum wires, where the configuration space is a graph, 
i.e. a strictly one-dimensional object and the Hamiltonian is minus the Laplacian with arbitrary 
boundary conditions at the vertices of the graph and which makes it a self-adjoint operator. The 
graph need not to be planar and may be bent when realized as a subset of the 3-dimensional 
Euclidean space M.^. By now many explicit examples have been considered (see e.g. [|l^, |2^, 



29, 3C, 31 



32, Bl, 34, 35, 36, 37]) including also the Dirac operator with suitable boundary 



conditions [38]. Our approach gives a systematic discussion and covers in particular all these 



cases for the Laplace operator. In this article we will, however, not be concerned with the 
question, which of these boundary conditions could be physically reasonable. The physical 



relevance of different boundary conditions is discussed e.g. in [31 



The scattering theory for these operators exhibit a very rich structure (see e.g. [4^, 41, 42]) 
and by Landauer's theory [ p3| ] provides the background for understanding conductivity in meso- 
scopic systems. The on-shell S-matrix at energy E is a.nn x n matrix if the graph has n open 
ends, which we will show to be given in closed matrix form in terms of the boundary condi- 
tions and the lengths of the internal lines of the graph. We exhibit covariance and invariance 
properties and show in particular that the on-shell S-matrix is symmetric for all energies if the 
boundary conditions are real in a sense which we will make precise. The main result of this 
article is that this on-shell S-matrix is unitary, continuous in the energy and even real analytic 
except for an at most denumerable set of energies without finite accumulation points. This set 
is given in terms of the boundary conditions and the lengths of the internal lines. This result 
may be viewed as the quantum version of Kirchhoff 's rule. For explicit examples this has been 
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known (see e.g. [|^, ^^), but again our approach provides a unified treatment. Pliysically tliis 
unitarity is to be expected since there is a local Kirclihoff rule. In fact, the boundary conditions 
imply that the quantum probability currents of the components of any wave packet associated to 
the different lines entering any vertex add up to zero. Our discussion of the boundary conditions 
will be based on Green's theorem and will just reflect this observation. We will actually give 
three different proofs, each of which will be of interest in its own right. 

Finally there is a general duality transformation on the boundary conditions (turning Dirich- 
let into Neumann boundary conditions and vice versa) which combined with an energy depen- 
dent scale transformation on the lengths of the internal lines relates the high energy behaviour 
of one theory to the low energy behaviour of the transformed theory. 

This article is organized as follows. In Section 2 we discuss the simple case with one vertex 
only but with an arbitrary number of open lines ending there. This will allow us to present 
the main elements of our strategy, which is the general theory of self-adjoint extensions of 
symmetric operators and its relation to boundary conditions in the context of Laplace operators. 
This discussion uses some elementary facts about Hermitian symplectic forms. Although some 
results will be proven again for the general set-up in Section 3, for pedagogical reasons and 
because they are easier and more transparent in this simple case, we will also give proofs. 
In Section 3 we discuss the general case with the techniques and mostly with proofs, which 
extend those of Section 2. We start with a general algebraic formulation of boundary conditions 
involving a finite set of half lines and finite intervals but without any reference to local boundary 
conditions on a particular graph. At the end we show that any of these boundary conditions may 
be viewed as local boundary conditions on a suitable (maximal) graph. 

The connection between the theory of self-adjoint extensions of symmetric operators and 
Hermitian symplectic forms was brought to the attention of one of the authors (R.S.) by G. 
Segal in 1987. In a recent paper S.P. Novikov stated that he learned this from I.M. Gelfand 
back in 1971. Unfortunately we were unable to trace back the precise history of this connection. 
It seems that it was made by several researchers at different times (see e.g. [ |46| , £7, 48, 49]) but 
still was not analyzed systematically. In Sections 2 and 3 and in Appendix A we will try to fill 
this gap. 

In Section 4 we consider the question what happens if one decomposes a graph into two 
or more components by cutting some of its internal lines and replacing them by semi-infinite 
lines. One would like to compare the on-shell S -matrices obtained in this way with the original 
one. If the graph has two open ends and and its subgraphs are connected by exactly one line, the 
Aktosun factorization formula [50] (see also [51, 52, 53, 54, 55, 56, 57, 58, 59,^]) for potential 
scattering on the line easily carries over to this case. Such rules are reminiscent of the Cutkosky 
cutting rules [61] for one-particle reducible Feynman diagrams. Also such relations are well 
known in network theories and then the composition law for the S-matrices figures under the 
name star product [51, pp. 285-286] (we would like to thank M. Baake (Tubingen) for pointing 
out this reference), [^]. If the cutting involves more than one line, the situation becomes more 
complicated and leads to interesting phenomena related to the semiclassical Gutzwiller formula 
and the Selberg trace formula [ p^ ] (see also [17, 18, |^]). We provide a general composition 
rule for unitary matrices, which we will call a generalized star product and by which the on- 
shell S-matrix of an arbitrary graph (with local boundary conditions) can be factorized in terms 
of the S-matrices of its subgraphs. We expect that this general, highly nonlinear composition 
rule could also be of relevance in other contexts. Note that there is some similarity between our 
results and the recursive approach of [62]. Section 5 contains a summary and an outlook for 
possible applications and further investigations. 

When this article was already submitted for publication we have learned about the work [ p^ ] 
(some results of which were announced in a short note [^]) by S.P. Novikov, where a program 
similar to ours is carryied out for discrete (combinatorial) Laplacians on tree graphs. 
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2 The Quantum Wire with a Single Vertex 



In this section we will consider a quantum wire with n open ends and joined at a single vertex. 
This toy model will already exhibit most of the essential features of the general case and is also 
of interest in its own right. In particular the general strategy and the main techniques of our 
approach will be formulated in this section. Let the Hilbert space be given as 

n = eti^i = ei=iL\[o, oo)). 

Elements £ H will be written as (V'l, ^^2, V'n) and we will call ipj the component of in 
channel j. The scalar product in H is 

n 
i=l 

with the standard scalar product on L'^{[0, oo)) on the right hand side. We consider the symmet- 
ric operator A*^ on H, such that 



with domain of definition V{I!\^) being the set of all where each ipi, 1 < i < n together 
with their first and second generalized derivatives belong to L^(0, oo) (i.e. ipi G Ty^'^(0, oo), a 
Sobolev space) and which vanish at x = together with their first derivatives. It is clear that A" 
has defect indices (n, n), such that the set of all self-adjoint extensions can be parametrized (in 
a noncanonical way) by the unitary group U(n), which has real dimension (see e.g. [ ]64| , 65]). 

There is an alternative and equivalent description of all self-adjoint extensions in terms of 
symplectic theory and which goes as follows. Let V C H he the set of all ip such that each 
'^i, ^ < i ^ n belongs to VF^'^(0, oo) - and we will then say ijj G VF^'^. On V consider the 
following skew-Hermitian quadratic form given as 



V) = (A</.,V^) - (0, AV') = -^{^A) 

with the Laplace operator A = d^ /dx"^ considered as a differential expression. Obviously 
Q vanishes identically on D(A'^). Any self-adjoint extension of A*^ is now given in terms 
of a maximal isotropic subspace of V, i.e. a maximal (linear) subspace on which 17 vanishes 
identically. This notion corresponds to that of Lagrangean subspaces in the context of Euclidean 
symplectic forms (see e.g. []66|]). To find these maximal isotropic subspaces we perform an 
integration by parts (Green's theorem) and obtain with ' denoting the derivative and ~ complex 
conjugation 



i=l 



We rewrite this in the following form. Let []■.!)—>■ C^" be the surjective linear map which 
associates to ip and tjj' their boundary values at the origin: 

Here T denotes the transpose, so [i/j], ip{0) and ip'{0) are considered to be column vectors of 
length 2n and n respectively. The kernel of the map [ ] is obviously equal to D(A'^). Then we 
have 
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where ( , )^2n now denotes the scalar product on C^" and where the 2n x 2n matrix J is the 
canonical symplectic matrix on C^": 

Here and in what follows I is the unit matrix for the given context. Note that the Hermitian 
symplectic form u) differs from the Euclidean symplectic form on C^" 

To find all maximal isotropic subspaces in V with respect to Q. it therefore suffices to find 
all maximal isotropic subspaces in C^" with respect to uj and to take their preimage under the 
map [ ]. The set of all maximal isotropic subspaces corresponds to the Lagrangean Grassmann 
manifold in the context of Euclidean symplectic forms []66[]. Since J is non degenerate, such 
spaces all have complex dimension equal to n. This description is the local Kirchhoff rule 
referred to in the Introduction. Moreover with M.^ denoting the orthogonal complement (with 
respect to (•, •)c2n) of a space we have the 

Lemma 2.1. A linear subspace Ai ofC'^"' is maximal isotropic ijf Ai^ = Jjvl and ijf Ai^ is 
maximal isotropic. 

The proof is standard and follows easily from the definition and the fact that = —I and 
Jt = -J. We use this result in the following form. Let the linear subspace M. = M.{A, B) of 
C^" be given as the set of all [-0] in C^" satisfying 

^VW + S-^'W = 0, (2) 

where A and B are two n x n matrices. If the n x 2n matrix {A, B) has maximal rank equal 
to n then obviously M. has dimension equal to n and in this way we may describe all subspaces 
of dimension equal to n. Also the image of C^" under the map {A, B) is then all of C" because 
of the general result that for any linear map T from C^" into C" one always has dim Ker(r) + 
dim Ran(T) = 2n. Writing the adjoint of any (not necessarily square) matrix X as X'^ = X^ 
we claim 

Lemma 2.2. Let A and B be be two n x n matrices such that {A,B) has maximal rank. Then 
A4{A, B) is maximal isotropic ijf AB'^ is self-adjoint. 

The proof is easily obtained by writing the condition (^ in the form (<I>'^, [V'Dc^™ = 0, 1 < 
k < n, where is given as the A;-th column vector of the 2n x n matrix (^4, S)-^ = {A, B)'^ . 
Obviously they are linearly independent. Then by the previous lemma M.{A, B) is maximal 
isotropic iff the space spanned by the is maximal isotropic. This condition in turn amounts 
to the condition that (yl, B)J{A, B)"^ = 0, which means that AB'^ has to be self-adjoint. The 
converse is also obviously true. 

Example 2.1. (n > 3, n odd) Consider 

i;j{0) + c,(V';_i(0) + V';+i(0)) = Ofor 1 < j < n 

with a mod n convention. The resulting A is the identity matrix and AB'^ is self-adjoint iff 
all Cj are equal ( =c) and real. Then B is of the form Bj^ = c{5jj^2 k + k+2)- (^fi^ even the 
condition is that Cj = Cj+i and Cj = Cj^2for all j) 

Example 2.2. (n > 3) Consider 

^/jjiO) + V,+i(0) + c,(V;(0) - V;+i(0)) = Ofor l<j<n 

again with a mod n convention. The resulting A has maximal rank and AB'^ is self-adjoint 
iff now all Cj are equal and purely imaginary. 
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Example 2.3. (n = 2) (see e.g. / p4 p3| pq/J To relate this case to familiar examples we realize 
H as L^(M) = cxD, 0]) © -^^^([0, oo)) and write the boundary conditions as 

V^'(0+) ; \c d )\ i;'{0- 



Then AB^ is self-adjoint iff the matrix 



a h 
c d 



belongs to SL(2, M) and fi is real. Up to a set of measure zero in U(2) this gives all self-adjoint 
extensions. The interpretation of the parameters entering the boundary conditions can be found 



in [68]. The case a— l = d— l = h = ^, exp(2z/x) = 1 gives the familiar 6-potential of 
strength c at the origin. The case a — l = d — l = c = 0, exp(2i^) = 1 gives the so called 
6' -interaction of strength b (see e.g. 57 ^ 7 <^nd references therein). The case a = d~ 



b = c = 0, exp(2i//) = 1 gives the 6' -potential of strength —2(1 — a) /(I + a) / |74 p^ . In 
particular for the choice a — l = d— l = b = c = 0, exp(2i//) = 1, the corresponding free 
S-matrix (see below) is given by 



1 

1 



More generally for n lines f~ Mj with free propagation and with the appropriate labeling of the 
2n ends the 2n x 2n on-shell S-matrix takes the form, 



St%E) = ^l Ij. (3) 

We will make use of this observation in Section 4, where we will exploit the fact that this 
^2n'^ (-B) serves as a unit matrix with respect to the generalized star product. 

In what follows we will always assume that A and B define a maximal isotropic subspace 
M. [A, B), such that the resulting operator, which we denote by A(^, B), is self-adjoint. A core 
of this operator 'D{/S.{A, B)) is given as the preimage of M.{/S.{A, B)) under the map [ ]. Note 
that D(A'^) has codimension equal to n in any of these cores. Thus the quantum mechanical one 
particle Hamiltonians we will consider are of the form — A(A, B) for any boundary condition 
[A, B) defining a maximal isotropic subspace. 

The self-adjointness of AB'^ , i.e. the relation AB^ — BA^ = 0, will be the main Leitmotiv 
throughout this article, since combined with the maximal rank condition it encodes the self- 
adjointness of the operator A{A, B) and is the algebraic formulation of the local Kirchhoff rule. 
The proof of Lemma 2.2 combined with the previous lemma also shows that M.{A, B)^ = 
JM.{A, B) = M.{—B, A). Note that {A\B'^) does not necessarily define a maximal subspace 
if {A, B) does. As an example let H be self-adjoint and A invertible and set B = H{A^^)'^ . 
Then (yl, B) defines a maximal isotropic subspace since AB'^ = H, but A^B = A^H{A-^)K 
B^A = A~^HA and these two expressions differ if AA"^ and H do not commute. If ^ = 
such that B is invertible we have Neumann boundary conditions and if B = such that A is 
invertible we have Dirichlet boundary conditions. 

We will now calculate the on-shell S-matrix S{E) = Sa,b{E) for all energies E > 0. This 
will be an n X n matrix whose matrix elements are defined by the following relations. We look 
for plane wave solutions E), 1 < k < not the Schrodinger equation for —A{A, B) in the 
form 

^^(xi?) = |^^-^^^^^^^' (4) 
' ' |e-*^^ + 5fcfc(^)e*^^' fori = fc 
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and which satisfy the boundary conditions. Thus Skk{E) has the interpretation of being the 
reflection amplitude in channel k while Sjk{E) with j 7^ is the transmission amplitude from 
channel k into channel j, both for an incoming plane wave exp{—i\/^x) in channel k. This 
definition of the S -matrix differs from the standard one used in potential scattering theory [fT^, 



73p , where the equal transmission amplitudes build up the diagonal. In particular forn = 2 and 



general boundary conditions at the origin as described in Example 2.3 we have 

S{E) = 



R{E) Ti{E) 
T2{E) L{E) 

r- ic f a-i^/Eb- ^-d 2e^^ 



^ ' "7 V 2e-*^ -a-WEb-^ + d 



whereas 

'S'standard(-E') = S{E) 



1 \ / T^iE) R{E) 

1 J V L{E) T2{E) 



These S-matrices are unitarily equivalent iff A(A, B) is real (i.e. e^*'' = 1 such that A(yl, B) 
commutes with complex conjugation, see also below for a general discussion) and is invariant 
with respect to reflection at the origin. In the latter case Ti{E) = T2(E) and R{E) = L{E). 

We return to the Ansatz (^. After a short calculation the boundary conditions for the ijj^ 
take the form of a matrix equation for S{E) 

{A + iVEB)S{E) = -{A-i^B). (5) 

To solve for S{E) we will establish the following 

Lemma 2.3. For all E > both matrices {A + i\/EB) and {A — i\fEB) are invertible. 
Proof: Assume det(^ + iVEB) = 0. But then also 



det(^"f - iVEB^) = det{A + i^/EB) = 0, 
so there is x 7^ such that {A'^ — i^/EB'^)x = 0. In particular 

= ix, (A + iVEB)iA^ - i^B^)x) = (Atx, ^t^) + E{B^x, B^), 

where we have used the fact that AB^ is self-adjoint. But this implies that A^x = = 0- 
Hence {Acp + B(j)', x) = for all cj), (j)' E C". But as already remarked the range of the map 
{A, B) : C^" C" is all of C". Hence x = and we have arrived at a contradiction. This 
proves the lemma since the invertibility of ^ — i\fEB is proved in the same way. Also we have 
shown that AA^ + EBB^ is a strictly positive operator and hence an invertible operator on 
for all E > 0. If A (or B) is invertible there is an easier proof of the lemma. Indeed, assume 
there is x / such that {A + iy/EB)x = 0. Then we have A^^Bx = i/y/E ■ x- But A^^B is 
self-adjoint since AB'^ is and therefore has only real eigenvalues giving a contradiction. To sum 
up we have proved the first part of 

Theorem 2.1. For the above quantum wire with one vertex and boundary conditions given by 
the pair {A^ B) the on-shell S-matrix is given as 

Sa,b{E) = -{A + i^/EB)-\A-iy/EB) 

= -{A^ -i^B^){AA^ +EBB^)-\A-i^B) 

and is unitary and real analytic in E > 0. 
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To prove unitarity, we observe that S^E) = -{A^ + ^^/EB'f){A'< - i^fEB'^)~'^ and 
S{E)~^ = —{A — i\/^B)~^{A + i^/EB). Now it is easy to see that these expressions are 
equal using again the fact that AB'^ is self-adjoint. Note that unitarity follows from abstract 
reasoning. In fact the difference of the resolvents of the operators /S.{A = 0,5 = 1) (Neumann 
boundary conditions and with S-matrix equal to I) and /\{A,B) by Krein's formula (see e.g. 



[|64|]) is a finite rank operator since we are dealing with finite defect indices, so in particular this 
difference is trace class. Therefore by the Birman-Kato theory the entire S-matrix exists and is 
unitary. Also the S-matrix for Dirichlet boundary conditions {A = 1,B = 0) is —I. The second 
relation in (^ can be used to discuss the bound state problem for — A(yl, B), which, however 
we will not do in this article. 

Relation (^ is a remarkable matrix analogue of the representation in potential scattering 
theory of the on-shell S-matrix at given angular momentum I as a quotient of Jost functions, i.e. 



Si{k = V E) = fi{k)/fi{—k) (see [74, 75]). As in potential scattering theory the scattering 
matrix Sa,b{E) as a function of \/E can be analytically continued to a meromorphic function 
in the whole complex plane. In fact, by the self-adjointness of IS.{A,B) it is analytic in the 
physical energy sheet Im \fE > except for poles on the positive imaginary axis corresponding 
to bound states and may have additional poles (i.e. resonances) in the unphysical energy sheet 
Im ^/E <0. 

If C is any invertible nxn matrix then obviously {CA, CB) defines the same boundary con- 
ditions as {A, B) such that M{CA, CB) = M{A, B) and A(C^, CB) = A(^, B) and this 
is reflected by the fact that Sca,cb{E) = Sa,b{E). Conversely, if M{A, B) = M{A', B') 
then there is an invertible C such that A = CA',B = CB'. This follows easily from the 
defining relation (^). We want to use this observation to show that the on-shell S-matrix 
uniquely fixes the boundary conditions. In fact assume that Sa.b{E) = Sa',b'{E) or equiv- 
alently Sa,b{E)Sa',b'{E)^ = I holds for some E > 0. By a short calculation this holds iff 
A'B^ - B'A^ = 0, i.e. {A', B')J{A, S)t = 0. But by the proof of Lemma 2.2 this means that 
M{A, B)^ = M{A', B')^ so the two maximal isotropic subspaces M{A, B) and M{A' , B') 
are equal as was the claim. 

To fix the freedom in parametrizing a maximal isotropic subspace by the pair [A, B), con- 
sider the Lie group Q{2n), consisting of all 2n x 2n matrices g which preserve the Hermitian 
symplectic structure, i.e. g'^ Jg = J. We claim that this group is isomorphic to the classical 
group U(n, n) (see e.g. [f?^). Indeed iJ is Hermitian, (iJ)^ = I and Tr J = such that the 
only eigenvalues ±1 are of equal multiplicity. Therefore there is a unitary W such that 



-I 

Thus elements g in the group WQ{2n)W~^ satisfy 

^ / I \ / I 

such that WQ{2n)W~^ = U(n, n). The group U(n,n) and hence Q{2n) has real dimension 
4n^ and the latter acts transitively on the set of all maximal isotropic subspaces. In particular 
Ai {A, B) is the image of M.{A = B = 1) under the map given by the group element 



B^{AA^ + BB^)-i A^{AA^ + BB^)-i 
-A^{AA^ + BB^)-l B^{AA^ + BB^)-\ 



Let/C(2n) be the isotropy group oi M.{A = Q,B = I), i.e. the subgproup which leaves M.{A = 
0, i? = I) fixed. It is of real dimension 3n^. The set of all maximal isotropic subspaces is in 
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one-to-one correspondence with the right coset space Q{2n) /K,{2n) which has dimension n? as 
it should. Also this space may be shown to be compact. In Appendix A we will relate the present 
discussion of selfadjoint extensions of A'^ with von Neumann's theory of selfadjoint extensions. 
In this context it is worthwhile to note that the parametrization of self-adjoint extensions in 
terms of maximal isotropic spaces is much more convenient for the description of these graph 
Hamiltonians, rather than the standard von Neumann's parametrization. In particular the content 
of Appendix A is not needed for an understanding of the main material presented in this article. 

Now we establish some properties of these on-shell S-matrices. Although we shall prove 
corresponding results in the general case they are more transparent and easier to prove in this 
simple situation. In particular, if the boundary condition {A, B) is such that A is invertible one 
may choose C = A^^ and similarly C = B^^ if B is invertible. To determine Sa,b{E) for all 
E in these cases it therefore suffices to diagonalize the self-adjoint matrices A^^B and B~^A 
respectively. Thus if A~^B = VHV~^ with a unitary V and diagonal, self-adjoint H, then 
Sa,b{E) = -V{I + iVEH)-^{I - iVEH)V-^. 

Thus in Example 2A H and V are given as Hki = 6ki2ccos 27r(/ — l)/n and 



resulting in an on-shell S-matrix of the form 



^ l + 2icVEcos^{l-l) 



1= 



With the equivalence {A, B) ~ {CA, CB) for invertible C in mind, it follows easily that 
Sa,b{E) is diagonal for all E iff A and B are both diagonal (Robin boundary conditions). 

Let [/ be a unitary operator on C^. Then U defines a unitary operator on in a natural 
way via {Uip)i = J2j=i ^iji^j- ^ special case this covers the situation where f7 is a gauge 
transformation of the form tpj —>■ exp{ixj)'ipj with constant Xj- Also {AU,BU) defines a 
maximal isotropic subspace and we have A{AU, BU) = U^^ A{A, B)l/(. Correspondingly 
we have Sau.bu{E) = U^^Sa,b{E)U. Assume in particular that there is a unitary U and 
an invertible C such that CA = AU and CB = BU. Then the relations A{AU,BU) = 
U-^A{A,B)U = A{A,B) and Sa,b{E) = U~^Sa,b{E)U are valid. This has a special 
application. There is a natural unitary representation tt ^ U (vr) of the permutation group of n 
elements into the unitaries of C" and analogously a representation vr — > U{-k) into the unitaries 
of H.. For any permutation vr such that AU{t^) = CA and BU (vr) = CB holds for an invertible 
C = C(7r) one has Sa,b{E) = U{'ir)^^ Sa,b{E)U{tt). In Examples 2.1 and 2.2 one may take 
TT to be the cyclic permutation and C(7r) = U{it). Consequently the on-shell S-matrix for the 
Example 2.1 given by satisfies Sj-^-i k+i{E) = Sjk{E) mod n for all I (see also e.g. [ p^ ] 
for other examples). 

Next we claim that if A and B are both real such that A'^ = A^ and = then Sa,b{E) 
equals its transpose. In particular the on-shell S-matrix (^ for the Example 2.1 is symmetric. 
This is in analogy to potential scattering on the line, where the Hamiltonian is also real (i.e. 
commutes with complex conjugation). There the transmission amplitude for the incoming plane 
wave from the left equals the transmission amplitude for the incoming plane wave from the 



right (see e.g. [ |72| , 73]). In fact one has a general result. For given {A,B) defining the self- 
adjoint operator A{A,B), {A,B) also defines a self-adjoint operator A{A,B), since {A,B) 
has maximal rank and AB'^ is also self-adjoint. It is easy to see that A{A,B) has domain 
V{A{A,B)) = {V; I G V{A{A,B)} and A{A,B)i) = A{A,B)i>. In particular A{A, B) 
is a real operator if the matrices A and B are real. Thus the Laplace operators obtained from 
Example 2. 1 are real while those obtained from Example 2.2 are not real. In Example 2.3 {A, B) 
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is real iff exp(2i//) = 1. More generally we will say that the boundary conditions {A, B) are real 
if there is an invertible map C such that CA and CB are real. Thus Robin boundary conditions 
are real. For invertible A (or B) a necessary and sufficient condition is that the self-adjoint 
matrix A^^B (or B~^A) is real. We have the 

Corollary 2.1. The on-shell S-matrices satisfy the following relation 

S^^siEf = Sa,b{E). (8) 

For the proof observe that S^ ^iEV = -(^^ - iVEB'f){A^ + i^/EB'^)-^ and the claim 
again easily follows from the fact that AB'^ is self-adjoint. 

As a next observation we want to exploit the fact that A and B play an almost symmetric 
role. First we recall that if {A,B) defines a maximal isotropic subspace then {—B,A) also 
defines a maximal isotropic subspace, which is just the orthogonal complement. In particular 
Dirichlet boundary conditions turn into Neumann boundary conditions and vice versa under 
this correspondence. Although the Laplace operators /S.{A,B) and /S.{—B,A) are not directly 
related, there is a relation for their on-shell S-matrices. In fact from ^ we immediately obtain 

Corollary 2.2. The on-shell S-matrices for the operators — A(j4, B) and —A{—B, A) are re- 
lated by 

S-ba{E) = -Sa,b{E-^) (9) 

for all E>0. 



Adapting the notation from string theory (see e.g. ||77|]) we call the map 9 : {A, B) ^ 
{—B^A) a duality transformation and (^ a duality relation. On the level of Laplace operators 
and on-shell S-matrices 6 is obviously an involution. Since M.{—B, A) = M.{A, B)-^ there are 
no self dual boundary conditions. 

We conclude this section by providing a necessary and sufficient condition for Sa,b{E) to be 
independent of E, i.e. to be a constant matrix. In fact this will occur if {A, B) is such that none 
of the linear conditions involve both and V''(0)- This means that the boundary conditions 
are scale invariant, i.e. invariant under the variable transformation x X^^x, x, A G M+. The 
algebraic formulation is given by 

Corollary 2.3. Assume the boundary condition is such that to any A > there is an invertible 
C(A) with C{X)A = A andC{X)B = XB. Then both Sa,b{E) and S-b,a{E) are independent 
of E. The only eigenvalues of Sa,b{E) are +1 and —1 with multiplicities equal to Rank B 
and Rank A. Also AB'^ = and both Sa,b{E) and S-b,a{E) are Hermitian such that by 
unitarity they are also involutive maps. If in addition the boundary conditions are real, then 
Sa,b{E) and S-b,a{E) are also real. Conversely if Sa,b{E) is constant then there is C(A) 
with the above property. 

Proof: The first part is obvious. To prove the second part recall (^). Since Sa,b{E) is 
constant this implies Sa,b{E)'^ B'^ = B^ and Sa,b{E)'^ A^ = —A"^. Thus the column vectors 
of B'^ and A^ are eigenvectors of SA.BiE^ = SA,BiE)~^ and therefore of Sa,b{E) with 
eigenvalues +1 and —1 respectively. They span subspaces of dimensions equal to Rank B^ = 
Rank B and Rank A"^ = Rank A. To see that these vectors combined span all of C" assume 
there is a vector if; orthogonal to both these spaces. But this means that Atjj = Bip = and this is 
possible only if = since A — i\/^B is invertible. Also eigenvectors for different eigenvalues 
are orthogonal which means AB"^ = 0. The hermiticity of Sa,b{E) and S-b,a{E) therefore 
follows easily from (^. The reality of Sa,b{E) and S-b,a{E) if the boundary conditions are 
real follows from Corollary 2.1. As for the converse we first observe that by the above arguments 



10 



a constant Sa^b[E) implies the previous properties concerning eigenvalues and eigenvectors. 
Let therefore C/ be a unitary map which diagonalizes Sa,b{E), i.e. 



USaAE)U-^ 



-I 
I 



holds with an obvious block notation. Also trivially 

USa,b{E)U-^UA^ = -UA\ USa,b{E)U-^UB^ = UBK 
Therefore A and B are necessarily of the form 



^21 



U, B 



b\ 
6'. 



12 



22 



u, 



again in the same block notation. Since {A, B) and hence {AU, BU) has maximal rank, the 
matrix 



*21 "22 



h' 



12 



is invertible. Let D denote its inverse, such that 



DA 



I 




U, DB 



Then 



C(A) = D 



-1 



I 
AI 




I 



D 



U. 



does the job, concluding the proof of the corollary. 

The following example will be reconsidered in Section 4. 

Example 2.4. (n = 3) Let the boundary conditions be given as 

V-lW =^-2(0) =^3(0), 

V'i(0) + V'2(0) + V'3(0) = 0, 



I.e. 



A 



1 


-1 

















1 




, B = 
















/ 




; 


1 


1 



, C(A) 




Thus A has rank 2, B rank 1 and AB'^ = 0. The resulting on-shell S-matrix is given as 

/ _1 2 

3 3 

SiE)= i ■ 

2 2 _i ; 

\ 3 3 3 / 



3 » 

1 2 
"3 3 



for all E. Since it is Hermitian and unitary the eigenvalues are indeed ±1 with the desired 
multiplicities since Ti S{E) = — 1. 
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3 Arbitrary Finite Quantum Wires 



In this section we will discuss the general case employing the methods used in the previous 
section. Let £ and T be finite sets with n and m elements respectively and ordered in an arbitrary 
but fixed way. £ labels the external lines and T labels the internal lines, i.e. we consider a graph 
with m internal (of finite length) and n external lines. Unless stated otherwise n / 0, since we 
will mainly focus on scattering theory. The discussion in the previous section already covered 
the case m = so we will also assume m 7^ 0. To each e G <S we associate the infinite interval 
[0, 00) and to each i £ Z the finite interval [0, Oj] where Cj > 0. The Hilbert space is now 
defined as 

'H = ©eefWe ©iGJ 'Hi = He © Ti-j, 

where He = ^^([0, 00)) for e G <S and Hi = ^^([0, flj]) for i £ T. He is called the exterior and 
Hj is called the interior component of H. Elements in H are written as 

Thus the previous case is the special case when I is empty. Let now A*^ be the Laplace operator 
with domain of definition D(A'^) being given as the set of all t/j with ipe, e G £ belongin to 
M^^'^(0, 00) and vanishing at x = together with their first derivatives while i/ji belong to 
W^^'^(0, Oj) for all i G X and vanish at the ends of the interval together their first derivatives. 
Obviously A*^ has defect indices {n + 2m, n + 2m). To find all self-adjoint extensions let now 
T> be the set of all tp with ipe,e£ £ belonging to VF2'2(0, 00) while the ijji G 14^2'2(0, Cj), i G I. 
Also the skew-Hermitian quadratic form on I? is now defined as 

0(0,^) = (A0,V) - (c/>,AV^) = -fliM^. 

Let [ ] : D ^ £<2(n+2m) ^jjg surjective linear map which associates to each the element [^] 
given as 

(V'e(0)ee^,V'i(0)iex,V'i(ai)iej)'^ \-( t\ nm 

again viewed as a column vector with the ordering given by the ordering of £ and T. Obviously 
P(A'') is the kernel of the map [ ]. By partial integration we again obtain 

where now J is the canonical symplectic form on ([;;2(n+2m) ^j^g same form as in ([T]) and 
( , )c2{n+2m) Is now thc canonical scalar product on ([;;2{"+2m) jj^g formulation of the boundary 
condition is as in the previous section. Let now A and i3 be (n + 2m) x (n + 2m) matrices and 
let M{A, B) be the linear space of all [V'] in c2("+2™) such that 

A'il^ + B'4/ = Q (11) 

Then B) has dimension n + 2m iff the (n + 2m) x 2(n + 2m) matrix {A, B) has maximal 

rank equal to n + 2m. If in addition AB'^ is self-adjoint then M.{A, B) is maximal isotropic. 
The resulting self-adjoint operator will again be denoted by A(yl, i?). Again in what follows 
we will always assume that the boundary conditions {A, B) have these properties. 

Similar to the discussion in the case of single vertex graph (see Appendix A) we can re- 
late von Neumann's parametrization of self-adjoint extensions of A*^ to the matrices A and B. 
Corresponding details can easily be worked out and therefore we omit them. 
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To determine the resulting on-shell S-matrix, we now look for plane wave solutions Tp''{-,E), 
k ^ £ of the stationary Schrodinger equation for — A(^, B) at energy E > which satisfy the 
boundary conditions {A, B) and which are of the following form and which generahze 

S^k{E)e'^^'' for jGSJ^k 

^-iVEx ^ Skk{E)e'^'' for j e£,j = k (12) 



ip'y{x,E) 

^ ajk{E)e'^'' + (3jk{E)e-'^'^ for j G I. 



The aim is thus to determine the n x n matrix S{E) = Sa,b{E) and the m x n matrices 
a{E) = aA,B{E) and P{E) = [5a,b{E). The physical interpretation of the matrix elements 
Sjk{E) is as before and the matrix elements of a{E) and P{E) are 'interior' amplitudes. It is 
advisable to write matrices like A and i? in a 3 x 3 block form, i.e. 



^11 


A12 


Ai3 


A21 


A22 


A23 


A31 


A32 


A33 



where An is an n x n matrix, A12 and A13 are n x m matrices etc. Thus for example the 
matrices An, An, Bii,Bii,i = 2, 3 describe the coupling of the exterior to the interior. 
Correspondingly ip and ip' are written as (compare (p^) 

4(0) 

± = I ^x(O) I , 3^' = I V'x(O) 




where a = (oi, ...am)- Also we introduce the diagonal m x m matrices exp(ibii/Ea:) by 

exp(±i\/^a)jfc = Sjue^'"^"^ for j, A; G T. 

The equations for the matrices S{E),a{E) and f3{E) now take the following form: 

SiE)+I \ / S{E)-I \ 

a{E) + P{E) + iVEB a{E) - (3{E) = 0, 

e'^^a{E) + e-^^^(3{E) J \ -e'^^a{E) + e-'^^(3{E) J 

which is a (n+2m) xn matrix equation with matrix multiplication between (n+2m) x (n+2m) 
and (n+2m) x n matrices. We rewrite this equation as an inhomogeneous equation generalizing 
equation (^, 

/ S{E) \ I I 

Za,b{E)\ a{E) \ =-{A-iVEB)\ | (13) 

V m J \ 

where 



with 



Za,b{E) = AX{E) + iVEBY{E) 



X{E) 





13 



and 



YiE) 



-e'^^ e 




If det 7^ the scattering matrix S{E) = Sa,b{E) as well as the m x n matrices 



a{E) and P{E) can be uniquely determined by solving the equation ([13[), 




-ZAMEyH^-i^^B) I 





(14) 



We recall that by the Birman-Kato theory Sa.b{E) is defined and unitary for almost all > 
because /S.{A,B) is a finite rank perturbation of A(A = 0,i3 = I). We denote by T,a,b = 
{E > I det Za,b{E) = 0} the set of exceptional points for which Za,b{E) is not invertible. 
Now we prove 

Theorem 3.1. For any boundary condition (A, B) the set Tia,b equals the set (ta,b of all pos- 
itive eigenvalues of — A(yl, B). This set is discrete and has no finite accumulation points in 
K+. 

Proof: First we prove that for every E G Sa.b all elements of Ker Za,b{E) have the form 
(0, 5?^, with S and [3 being (column) vectors in C". Let us suppose the converse, i.e. that 
there is a column vector (S^, 6F,0^)'^ G £n+2m ^^^j^ ^^^^ 



Za,b{E) I a 



(15) 



or equivalently 



{A + i^/EB) a \ +{A- i\/EB) 



( 



As in the proofs of Lemma 2.3 and Theorem 2.1 A + i\/ EB and A — i\/ EB are invertible and 
{A + iVEB)~'^{A - iVEB) is unitary such that 




-{A + iVEBy'^iA - iy/EB) 



( 

V 



Since unitary transformations preserve the Euclidean norm we get 



Nllc" + IIq^IIc'" + ll/^llc" 



Iq^IIc'" + ll/^llc™) 



such that s = 0. 



Now for arbitrary (0, a^,(^)'^ G Ker Za,b{E) we consider 

, . ^ f for j G £, 

^^^^ \ aje'^'' +I3je-''^'-' for j G T. 
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Obviously ijj{x) is an eigenfunction of —A{A, B). Thus we have proved that Ti^^b ^ ca,b- We 
note that this inclusion is nontrivial, since a priori we cannot exclude real energy resonances. 

Conversely, let E G (Ja,b- Observing that positive energy eigenfunctions must have support 
on the internal lines of the graph and repeating the arguments, which led to equation (|^, we see 
that there exists a nonzero vector (0, a^, € £n+2m ^^^j^ ^^mi (15) is satisfied, and hence 
<7A,B C We notice that det Za.b{E) is an entire function of \fE in the complex plane C. 

Also det Za.b{E) does not vanish identically since by the preceding arguments det Za,b{E) = 
for £^ G C \ M imphes that E is an eigenvalue of —/S.{A,B), which in turn contradicts the 
self-adjointness. Thus all real zeroes of det Za,b{E) are isolated. This concludes the proof of 
the theorem. 

Remark 3.1. The case n = with no associated S-matrix is of interest in its own right. Then 
( |T3| ) takes the form of a homogeneous equation 

Za,b{E) { ! 1=0, (16) 



I3{E) 

with a{E) and (3{E) being column vectors in C". It has solutions iff E €z M is such that 



det Za,b{E) = 0. By what has been said so far it is clear that the solutions of (16) give the 
eigenvalues and the eigenfunctions of a quantum wire without open ends. As a possible example 
one might consider the graph associated to the semiclassical description of the fullerene and 



choose appropriate boundary conditions at the vertices (see e.g. [78]). 



Before we proceed further with the study of the equation ( |13| ) we consider several examples. 
Starting with the case, where everything works well without any singularities, we consider the 
following 

Example 3.1. (n = m = 1) We choose the following realization of the Hilbert space 7i: 

n = L^([0,oo)) = L2([l,cx)))eL2([0,l]). 
We take Robin boundary conditions at the origin and a 6 potential of strength cat x = 1, i.e. 

sin (f V'(O) + cos (f ip'{0) = 

V^(i+)-V(i-) = 
^^'(1+) - V'(i-) -cV^(i) = . 

S{E), a{E) and P{E) are now functions and a straightforward calculation gives the following 
result. Let Sji{E) denote the on-shell S-matrix for the Robin boundary condition alone, i.e. 



Then 



Sr{E) = e2^^«(^) = sin(^-i^cosv9 ^ 
sm.(p + i^/E cos if 



(2^v^+c)e'(^+^'^-R(^»+ce-'^ -2i^ 



a{E) 



(2iv^-c)e-'^-ce'(v^+2*-R{£;)) 



(2jv^-c)e-'v^+ce'(v^+2'5fl(E)) 
2tV£e-'^ 



(2iv^-c)e-^v^+ce'(v^+2'5fl(-E)) ' 

In particular these quantities are finite for all E > 0. 
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3 . 



Figure 1: The graph from Example 3.2. The arrows show the positive direction for every 
segment. 



Example 3.2. Consider the graph depicted in Fig. with £ = {1,2}, X = {3, 4}, and with 
03 = 04 = a, i.e. the internal lines have equal length. The arrows on the graph segments show 
the positive directions. We pose the following real boundary conditions, obviously defining a 
self-adjoint operator, 

^^2(0) = ^3(0) = ^"4(0), 

V''i(o) + V'^(o)+vl(o) = 0, 
^2(0) - ^-3(0) - ^» = 0. 

Straightforward calculations yields 

detZ(S) = (10 - e^^^'^ _ 9e-2iv^«)^; 

2 2 

such that T,A,B = { ^Ji 1 ^ N}. The corresponding wave functions have the form ipi = 
tjj2 = 0, ip3{x) = —^^4(2;) = sin(^^). Note that these eigenfunctions for the embedded 
eigenvalues have compact support. This is in contrast to standard Schrddinger operators of 
the form —A.+V, where bound state eigenfunctions can not have compact support due to the 
ellipticity of —IS.. The on-shell S-matrix can also be easily calculated, giving 



S{E) 



1 



g^iVEa 3(e2iv^" - 1) 



(17) 



Equation (p^ is solvable for all E G T,a,b cind defines S{E) uniquely. S{E) is well behaved 
for all E > and there is no infiuence of the bound states, except that the (equal) refiection 
amplitudes vanish for E S '^a,b- This is in contrast to Schrddinger operators of the form 
— A + \P with appropriate real A and with P being the orthogonal projector onto any given 
one-dimensional linear subspace of the Hilbert space /|7^, or to Schrddinger operators with 
Wigner - von Neumann type potentials The scattering matrix (17) has poles only in the 

unphysical sheet (resonances). 



— 



.log 9 
2a 
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The relation 



holds as to be expected from the boundary conditions. 
Let us observe that 

det X{E) = det Y{E) = JJ(-2i sin \/^aj). 

Therefore if E G = UjexSa(j) = Ujex{-E > | sinVEaj = 0} then Za,b{E) will 
not be invertible for {A, B) defining Diriclilet or Neumann boundary conditions, since then 
det Za,b{E) = det X{E). In particular in these two cases the exterior and the interior 
decouple such that then we have 

A(yl, B) = Af (A, B) e Ai(yl, B) 

Here A£-(^, 5) for both {A = \B = and = 0, i? = I) have an absolutely continuous 
spectrum and Ai(j4, B) has a purely discrete spectrum which on the set > equals Eq. 
This means that we have eigenvalues embedded in the continuum. Now the equation (pj|) for 
Dirichlet or Neumann boundary conditions has a unique solution S{E) = —I and S{E) = I, 
respectively, and a{E) = P{E) = whenever E is not in Eq. If E is in Ea, then S{E) 
is still of this form but a{E) and f3{E) are nonunique and of the form a{E) = (3{E) and 
a{E) = —f3{E) for Dirichlet and Neumann boundary conditions, respectively, with ajk being 
arbitrary whenever E £ ^a{j) and zero otherwise. 

A similar result is vahd for arbitrary boundary conditions: 



Theorem 3.2. For any boundary condition {A,B) and all E G T,a,b th^ equation (13) is 
solvable and determines Sa.b{E) uniquely. 

Proof: Let ijj be an arbitrary element of Ker Z{E)\ i.e. ip solves the adjoint homogeneous 
equation corresponding to ([l3|), 

{X{E)'^A^ - i\/EY {E)"^ B^)i) = 0. (18) 

Suppose first that i? = 0, such that A has full rank (Dirichlet boundary conditions). Moreover 
T,A,B = Sq. Any solution of (18) may be written d& = A"^ x with x being a solution 



of X{E)'^x = 0- X is necessarily of the form x = (0, a^,/?^)"^ for some S,/? G C". Let 
(pi G C""*"^™, i = 1, . . . , n be such that {(pi)^ = ^ik- Then 

{^,ApP) = {A^xP,Pi) = {x,Pi) = Q 

for alH = 1, . . . ,n. We now apply the Fredholm alternative as follows. First we note that 
Ker Z{E)'^ = (Ran Z{E))-^, the orthogonal complement of the range of Z{E). Therefore the 
last relation states that all column vectors of 

{A - i^B) 



are in the range of Z{E). Thus ( [T3| ) has a solution and since all elements of Ker Z{E) are of 
the form (0, , the on-shell S-matrix Sa,b{E) is determined uniquely. 
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Now we suppose that B ^ 0. First we consider the case E e 'Ea,b and E ^ E<j. From (|l|) 
it follows that 

= iVEX{E)^~W(E)'^B^i), 

and thus 

(-0, BA^) = iVE{B^i), X{E)^~W{E)^B^i;). (19) 
Since BA^ is self-adjoint the l.h.s. of (19) is real. To analyze the r.h.s. we note that 



X{Ey Y{E) 



10 





/ 

cotan(\/^a) 







1 







cotan(\/^a) 



sin(v^a;) 

where the first term is self-adjoint and the second skew-self-adjoint. Therefore 

Re{B''tp,X{E)^~W{E)''B''^p) = i B'fijA 1 







Since the l.h.s. of (|^) is real, it follows that 



and thus 



I 
B^i/j, I I B^ 




= 



(20) 



(21) 



(22) 



for alH = 1, . . . , n. Multiplying the equation ( |18[ ) by 0j from the left we obtain 

{(f>i,A^ij)-iVE{cl)i,B^) = 0. 
Thus from (22) it follows that A^il^) = 0. Therefore 

(ip, {A - i\fEB)(t)i) = 

for alH = 1, . . . ,n. Therefore we may again invoke the Fredholm alternative and obtain a 
unique on-shell S-matrix. Finally we turn to the case E G Sa.b H Sq. An important ingredient 
of the proof in this case is the Moore-Penrose generalized inverse (or pseudoinverse) (see e.g. 
||8l|]). Recall that for any (not necessary square) matrix M its generalized inverse M* is uniquely 
defined by the Penrose equations 

MM*M = M, M*MM* = M*, 
{M*M)^ = M*M, (MM*)^ = MM*. 



Also one has 



(M*)t 
Ran M* 
Ker M" 



(M^r, 

RanM^ 
KerM'f, 
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and MM* = -PRanA/> M*M = -P^anA/t' where P-j-i denotes the orthogonal projector onto the 
subspace H. Moreover, 0* = 0. If M is a square matrix of full rank M* = M^^. However, 
the product formula for inverse matrices (M1M2)* = M* does not hold in general. The 
pseudoinverse of any diagonal matrix A with det A = is given by 



A* 



v 



\ 



0/ V 



/ 



Using the fact that {Q^ MQY = Q^M*Q for any unitary Q one can easily calculate the gener- 
alized inverse by means of the formulas 

M* = {M^MyU^ = M\MM^y. 

With these preparatory remarks we return to the proof. We will say that x ^ KerX(£')1^ 
is a basis element of Ker X{E)'^ if X = (0) oF ,1^^)^ with Sj = /3i = for alH = 1, . . . , m 
except for some k = 1, . . . ,m and 3^ / 0, / 0. For any basis element x either 



X{E)x = and Y{E)x = cyx 



(23) 



or 



Y{E)x = and X{E)x = cxX (24) 
with some cx , cy 7^ 0. The proof is elementary and is left to the reader. Taking the scalar 



product of ( 18) with x we obtain 

{X{E)x,A^i^) = iVE{Y{E)x,B^). 

In the case ( ^3| ) we have (x, B'^iI))=0, and in the case (|2^ {x, A'^^|J) = 0. 
Multiplying equation ([l8|) by X{E)*'^ we obtain 



Pi 



RanX{E) 



A^ip = i^/EX{Ey'^Y{E)^B^il:, 



and thus 



{^,BP^,,x{E)AH) = i^{B^i^,X{Ef^Y{EyBH)- 



(25) 



The matrix X{EY'^Y{E)'^ has the form of the r.h.s. of (2C), where the singular entries must be 
replaced by zeroes. Thus 

I 

Re{B^,X{Ey^Y{EYB^) = | ^t^, | | 5^-0 | . 



Using the relation P^^^ x(E) = I - ^Ran x{E)^ = I - Pkct x{Ey we may rewrite the l.h.s. of 
&) as 
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With Xi (1 < ^ < dim Ker X{E)^) being the basis elements of Kev X{E)^ we have 

dimKerX(£;)t 

1=1 



By the discussion above all terms in this sum are zero. Thus we again obtain (21). The rest 
of the proof is as in the preceding cases. Theorem 3.2 in particular says that the presence of 
bound states does not spoil the existence and uniqueness of the on-sheU S-matrix. This is to be 
expected since bound states do not participate in the scattering. The same is true for resonances 
due to their finite lifetime. 

We are now prepared to formulate the main result of this article 

Theorem 3.3. For any boundary condition {A, B) defining the self-adjoint operator — A(A, B) 
the resulting on-shell S-matrix Sa,b{E) is unitary for all E G ]R-|_. 

Proof: The proof is quite elementary. By arguments used in the proof of Theorem 3.1 we 
obtain that any solution of (pj|) satisfies 

SiE) \ / I 

a{E) \ = -{A + iVEB)-\A-iVEB) i f3{E) 
e-i^iL(3[E) I V e'^^a{E) 

Now {A + iy/EB)~^{A — i^/EB) is unitary. Multiplying each side with its adjoint from the 
left we therefore obtain 

S{E)'^S{E) + a{E)^a{E) + p{E)^ f3{E) = I + P{E)^ p{E) + a{E)^a{E) 

which gives S{E)^ S{E) = I and which is unitarity. A second 'analytic' proof is given in 
Appendix B (see also Section 4 for a third proof based on the generalized star product). 

We now discuss some properties of Sa,b{E). Obviously Za,b{E) can be analytically 
continued to the complex \/^-plane. By the self-adjointness of A{A,B) the determinant 
det Za,b{E) cannot have zeroes for E with Im ^fE > except for those on the positive imag- 
inary semiaxis corresponding to the bound states. It may have additional zeroes (resonances) in 
the unphysical energy sheet Im ^/E < 0. By means of (Q) the scattering matrix Sa,b{E) can 
be analytically continued to the whole complex plane as a meromorphic function of \/^- In fact 
it is a rational function in \/^, exp(i\/^aj), and ex.p{—i\/^aj). Therefore it is real analytic 
for all E G M+ \ T,a,b- Since Sa,b{E) for E E Tia,b is well defined and unitary, it is also 
continuous for all E £W. 

As for the low and high energy behaviour we have the following obvious property. If A is 
invertible then lim^jo Sa,b{E) = —I and if B is invertible then limEioo Sa,b{E) = I. For 
arbitrary {A, B) the corresponding relations do not hold in general as may be seen from looking 
at the Dirichlet {A = 1,B = {)) and Neumann {A = B = 1) boundary conditions. 

Let C be an invertible map on C"+2™ such that A(C^, CB) = A(A, B). Correspondingly 
we have Sca,cb{E) = Sa,b{E) as it should be, since Zca,cb{E) = CZa,b{E). Furthermore 
let [/ be a unitary map on C". Then U induces a unitary map U on by 

U = 

and a unitary U onH via 




^ipj otherwise 
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such that A{AU,BU) = U-^/\{A,B)U. We recall that all spaces H,, (e G 8) are the 
space i^^([0, oo)), so the definition oiU makes sense. Also Zj^^ bu^^) ~ ^a,b{E)U since U 
commutes with X{E) and Y{E), such that 



'^aO bu- ^^^^ observe that 



/ s 



u 



AU,BU 
^AU,BU 
V ^AU,BU 



(E) 
(E) 
(E) 



I US 



AU,BU 



AU,BU 
AU,BU 



(E) \ 

(E) 
(E) 



and 



{AU - iVEBU) 




{A - iVEB 




From this we immediately obtain first for E outside T,a,b and then by continuity for all > 
CoroDary 3.1. The following covariance properties hold for all E > 



^AU,BU 
^AU,BU 
^AU,BU 



{E) = U-^Sa,b{E)U 
{E) = aA,B{E)U 
{E) = I3a,b{E)U. 



(26) 



In particular if U is such that there exists an invertible C 
CB = BU then Sa,b{E) = U-^Sa,b{E)U for all E > 0. 



C{U) with CA = AU and 



We have the following special appUcation. There is a canonical representation t: U{7r) of 
the permutation group of n elements into the unitaries of C". If there is a tt and an invertible C = 
C(7r) such that CA = AU{Tr) and CB = BU{tt) then Sa,b{E) = U-^{-k)Sa,b{E)U{-k) for 
all E > {). The on-shell S-matrix in Example 3.2 is obviously invariant under the permutation 

1 •i-^ 2 of the two external legs. 

Remark 3.2. This discussion may be extended to the case that some of the interval lengths ai 
are equal resulting in more general covariance and possibly invariance properties described by 
some U, U and U such that Uu = U and Un = tin = Ofor i = 2,3. We leave out the details, 
which may be worked out easily. 

Corollary 2. 1 has the following generalization 

Corollary 3.2. For all boundary conditions {A,B) and all E > the following relation holds 



S,^{Ef = SA,B{E)- 



(27) 



In particular for real boundary conditions the transmission coefficient from channel j to channel 
k, k j {j,k & £) equals the transmission coefficient from channel k to channel j for all 
E>0. 

Proof: We start with two remarks. First the relation 



{A - iVeB) I ° j = ~ i^BY' 




holds for any X' of the form 



X' 



I 


















^'2,2 


^33 
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and similarly for Y'. Secondly the matrix 



aA,B{E) 
Pa,b{E) J 

constitutes the first n columns of the matrix -Za,b{E)~''-{AX' - i^/EBY'B) with X' and Y' 
arbitrary as above. Assume for a moment that E is not in S^. To prove the corollary for such E 
it therefore suffices to show that 

{ZA^s{Er\AY{E)-^^ - iVEBX{E)-^^)f 
= ZAMErHAYiE)-'^ - ^^/EBX{E)-'^), 

which is equivalent to 

{AX{E) + iVEBY{E)){Y{E)~^A^ - iVEX{E)'^B^) 
= {AY{E)-^'^ -i^/EBX{E)-^^){X{EfA^ + i^/EY{EfB^). 

But this relation follows from the self-adjointness of AB^ and the observation that X{E)Y {E)^^ 
is symmetric, such that the relations X{E)Y{E)-^ = Y{E)-^'^ X{EY and Y{E)X{E)-^ = 
X{E)~^'^Y{E)'^ hold. Finally by continuity we may drop the condition that E is not in S^, 
thus completing the proof. The on-shell S -matrix of Example 3.2 is obviously symmetric. 

To generalize the duality map of the previous section, we have to take the a dependence into 
account, so we write SA,B,a etc. The reason is that the length scales aj induce corresponding 
energy scales. Also the Hilbert spaces depend on a, 7i = TCa, so we will compare on-shell 
S -matrices related to theories in different Hilbert spaces. 

For given a let a{E) be given by ai{E) = Eai,i G X such that = Xa{E) and 

(E-^) = Ya{E) for all E > 0. Also set 

/ I \ 

T= I , 

V -I J 

such that rt = r, = I and TXa{E)T = Ya{E). We now define 9{A, B) = {-BT, AT). 
It is easy to see that 9{A, B) defines a maximal subspace. Also E~^ is not in '^0[A,B),a{E) if E 
is not in TiA,B,a and vice versa. 

This leads to following generalization of Corollary 2.2 

CoroUary 3.3. For all boundary conditions (A,B) and all E > the following identities hold 

S0(A,B),a{E){E~^) = — SA,B,aiE) 

(^0{A,B),aiE){E~^) = - aA,B,a{E) (28) 
f^e{A,B),a{E){E~^) = PA,B,a{E). 

The proof is easy by observing that 

Ze{A,B),aiE){E~'^) = --^ZA,B,a{E)T 

and 

~BT - -^AT = — ^-={A - iy/EB)T. 
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We conclude this section by giving a geometrical description of an arbitrary boundary con- 
dition as a local boundary condition at the vertices of a suitable graph. 

For given sets £ and I and a, we label the halfline [0, oo) associated to e G by Ig = 
[Oe,oOe) and the closed interval [0, a^] associated to i G X by Jj = [Oi,aj] (considering Oj as 
a generic variable there should be no confusion between the number Oj and the label Oj). By 
/ = Ueeele Ujgi li we denote the disjoint union. Let V = Uee£-{Oe} Ujgi {Oi,ai} C / be 
the set of 'endpoints' in /. Clearly the number of points in V equals | £ \ +2 \ I \= n + 2m. 
Consider a decomposition 

V = U^esV^ (29) 

of V into nonempty disjoint subsets with S being just an index set. We say that the points 
in V C / are equivalent (~) when they lie in the same V^. By identifying equivalent points 
in V C / we obtain a graph T, T = 1/ ~. In mathematical language T is a one-dimensional 
simplicial complex, which in particular is a topological space and noncompact if £ is nonempty. 
Obviously the vertices in T are in one-to-one correspondence with the elements ^ E H. Note 
that r need not be connected. Also there may be 'tadpoles', i.e. we allow that Oj and Cj for some 
i G X belong to a same set Vg. There is no restriction on the number of lines entering a vertex. 



In particular this number may equal 1 (so called dead end side branches | IC]), see Example 3.1. 
The graph need not be planar. 

Let {A, B} be the equivalence class of the boundary condition {A, B) with respect to the 
equivalence relation given as (A', B') ~ {A, B) iff there exists an invertible C such that A' = 
CA, B' = CB. By our previous discussion /S.{A,B) only depends on {A,B}. We say that 
{A, B} has a description as a local boundary condition on the graph F if the following holds. 
First observe that we may label the columns of A and B by the elements w in V (see (|l^)). With 
this convention there is supposed to exist (^4', B') G {A, B} with the following properties. To 
each k labeling the rows of A' and B' there is ^ = £,{k), such that A'^^ = = for all v not 
in V^. In other words the boundary condition labeled by k only involves the value of if) and its 
derivative at those points in V which belong to and this set is in one-to-one correspondence 
with a vertex in T. Of course if F is the unique graph consisting of one vertex only then this 
F does the job for any boundary condition {A,B}. However, one may convince oneself that 
for any given boundary condition {A,B} there is a unique maximal graph F = T{{A,B}) 
describing {A, B} as a local boundary condition, where maximal means that the number of 
vertices is maximal. 

Let us briefly indicate the proof. Arrange the 2(n+2m) columns of the (n+2m) x 2(n+2m) 
matrix {A, B) in such a way that the first n + 2m columns are linearly independent. Call this 
matrix X. Then there is an invertible matrix C such that the (n+2m) x (n+2m) matrix made of 
the first n + 2m columns of CX is the unit matrix. Now rearrange CXhy undoing the previous 



arrangement giving (CA, CB).The decomposition (29) and hence F({A, B}) may be read off 
the 'connectivity' of {CA,CB). The self-adjointness of AB^ may now be verified locally at 
each vertex (the local Kirchhoff rule), see Examples 3.1 and 3.2. Of course different boundary 
conditions may still give the same F in this way. Also in this sense the graphs associated to the 
discussion in Section 2 may actually consist of several disconnected parts, each with one vertex 
but without tadpoles. 

In particular this discussion shows that the results of this section cover all local boundary 
conditions on all graphs with arbitrary lengths in the interior and which describe Hamiltonians 
with free propagation away from the vertices. 
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4 The Generalized Star Product and Factorization of the S-matrix 



In this section we will define a new composition rule for unitary matrices not necessarily of 
equal rank. This composition rule will generalize the star product for unitary 2x2 matrices 
so we will call it a generalized star product. It will be associative and the resulting matrix will 
again be unitary. We will apply this new composition rule to obtain the on-shell S-matrix for 
a graph from the on-shell S-matrices at the same energy of two subgraphs obtained by cutting 
the graph along an arbitrary numbers of Unes. By iteration this will in particular allow us to 
obtain the on-shell S-matrix for an arbitrary graph from the on-shell S-matrices associated to its 
vertices (see Section 2), thus leading to a third proof of unitarity. 

Let V be any unitary p x p matrix {p > 0). The composition rule will depend on V and 
will be denoted by * v , such that for any unitary n' x n' matrix U' with n' > p and any unitary 
n" X n" matrix U" with n" > p, 2p < n' + n" and subject to a certain condition (see below) 
there will be a resulting unitary n x n mattix U = U' *v U" with n = n' + n" — 2p. This 
generaUzed star product may be viewed as an amalgamation of U' and U" and with V acting as 
an amalgam. To construct *y we write U' and U" in a 2 x 2-block form 

where C/22 ^^^^ are p x p mattices, U[i is an (n' — p) x (n' — p) matrix, U22 is an (n" — 
p) X {n" — p) matrix etc. The unitarity condition for U' then reads 



U''^ J J' 


^ '^2l'^21 


= I 




+ ^22^22 


= I 


'^ll'-^12 


+ ^21^22 


= 


TT''^ TJ' 

'^12<^ll 


+ ^^22^^21 


= 



and similarly for U" . 

Condition A : The pxp matrix yc/22^~^^n "^oes not have 1 as an eigenvalue. 

Note that by unitarity of U' , U" and V one has || VU^iV-'^Ul^ ||< 1. Strict inequahty 
holds whenever || C/22 II < 1 oi" II ^ii II < 1 and then Condition A is satisfied. In general if 
Condition A is satisfied it is easy to see that the following pxp matrices exist. 

K2 = (I - V-^U'^^VU^2)~^V~^ = ^"Hl - ^^n^^^22^"^)"^- 
An easy calculation estabhshes the following relations 

Ki = V + VU'22V-^U'{^Ki = V + VU'22K2U'(^V 

V + KiU!22V-^U'{^V, 

(31) 

K2 = V-'^ + V~^U'{^VU'22K2 = V~^ + V~'^U'{^K]U'22V-^ 

V-^ + K2U'{^VU!22V-^ . 



24 



Note that formally one has 



m=Q 

oo 

K2 = Y.^^''^nyu^2ry-'- 



-=o (32) 



m=0 



With these preparations the matrix U = U' *vU" is now defined as follows. Write i7 in a 2 x 2 
block form as 



U 



Uu Ui2 

U21 U22 



(33) 



where Un is an (n' — p) x (n' — p) matrix, U22 is an {n" — p) x {n" — p) matrix etc. These 
matrices are now defined as 

= U[^ + U[^K2U1^VU!,^, 

U22 = U^2 + U^iKiUl,^V-W^^, 

U12 = U[^K2U1^, 

U21 = U^iKiUl,^. 

In particular if n' = 2p then 

I =( I r i I 



Similarly if n" = 2p then 



10 \ V \ 



In this sense the matrices ( ? !^ 1 serve as units when V = 1 



I 

A straightforward but somewhat lengthy calculation presented in Appendix C gives 

Theorem 4.1. If Condition A is satisfied then the matrix U = U' *v U" is unitary. 

Analogously one may prove associativity. More precisely let U'" be a unitary n'" x n'" and 
V a unitary p' x p' matrix with p' < n", p' < n'". If p + ^i' < n', then 

U' *v {U" *v' U'") = {W *v U") *v' U'" 

holds whenever Condition A is satisfied for the compositions involved. 

We apply this to the on-shell S -matrices of quantum wires as follows. For the special case 
F = I we introduce the notation *p = *v=i- Let T' and T" be two graphs with n' and n" 
external fines labeled by S' and £", i.e. \£'\ = n', \£"\ = n" and an arbitrary number of internal 
lines. Furthermore at all vertices we have local boundary conditions giving Laplace operators 
A(r') on r' and A(r") on T" and on-shell S-matrices S'{E) and S"{E). Let now £^ and S^^ 
be subsets of £' and £" respectively having an equal number {= p > 0) of elements. Also let 
ipo : £q ^ £q be a one-to-one map. Finally to each k e SqV/c associate a number ak > 0. 
With these data we can now form a graph T by connecting the external fine k e £q with the line 
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ipo{k) E £q to form a line of length Cfc. In other words the intervals [0^, oofc) belonging to F' 
and the intervals [0(pg(fc), oo<^„(,fc)) belonging to T" are replaced by the the finite interval [Ofc, a/c] 
with Ofc being associated to the same vertex in T' as previously and being associated to the 
same vertex in T" as O^^f^k^ before in the sense of the discussion at the end of Section 3. Recall 
that the graphs need not be planar. Thus F has n = n' + n" — 2p external lines indexed by 
elements in {£' \£q) U {£" \ £q) and p internal lines indexed by elements in £q in addition to 
those of F' and F". There are no new vertices in addition to those of F' and F" so the boundary 
conditions on F' and F" define boundary conditions on F resulting in a Laplace operator A(F). 
The following formula relates the corresponding on-shell S-matrices S'{E), S"{E) and S{E). 
First let the indices of £q in £' come after the indices in £' \ £q (in an arbitrary but fixed order) 
(see (|30|)). Via the map ipo we may identify £q with £q so let these indices now come first in £", 
but again in the same order. Finally let the diagonal matrix V{a) be given as 



V{a) 



exp i\fEa 
I 



where exp(f\/^a) again is the diagonal p x p matrix given by the p lengths a^, ke£'Q. Then we 
claim that the relation 

S{E) = S\E) *p V{a)S"{E)V{a) (34) 

holds. The operators Ki and K2 involved now depend on E and are singular for £^ in a denu- 
merable set, namely when Condition A is violated. This follows by arguments similar to those 



made after Theorem 3.3. Thus these values have to be left out in (|34t). In the end one may then 
extend (^^ to these singular values of E by arguments analogous to those after Remark 3.1. If 
F is simply the disjoint union of F' and F", i.e. if no connections are made (corresponding to 
p = and n = n' + n"), then S{E) is just the tensor product of S'{E) and S"{E). In this 
sense the generalized star product is a generalization of the tensor product. Also by a previous 
discussion (see (|3|)) V^^S{E)V = S2n^{E) *v S{E) for any on-shell S-matrix with n open 
ends and any unitary nxn matrix V. Similarly S{E) *v S2n''{E) = VS{E)V-'^. Using (g) 
the on-shell S-matrix associated to any graph and its boundary conditions is obtained from the 
on-shell S-matrices associated to its subgraphs each having one vertex only. In fact, pick one 
vertex and choose all the internal lines connecting to all other vertices. This leads to two graphs 



and the rule (34) may be applied. Iterating this procedure L times, where L is the number of 



vertices, gives the desired result. 




Figure 2: Decomposition of the graph. 
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We will not prove the claim (^) here but only give formal and intuitive arguments which 
also apply in the physical context of the usual star product or Aktosun formula for potential 
scattering on the line and which are based on a rearrangement of the Bom series for the on-shell 
S -matrix. For partial results in potential scattering in higher dimensions when the separation of 



two (or more) potentials tends to infinity see \ p2l |83| ]. The complete proof of (g4j) will be given 
elsewhere [|^. 

For the sake of definiteness we consider the amplitudes S{E)j.i, k,l £ £' \ £q, which form 
Uu in this context, since we have U' = S'{E) etc. The other amplitudes may be discussed 
analogously. The first nontrivial contribution is S'{E)ki corresponding to the first term in the 



first relation of (33). In next order the incoming plane wave in channel I may cross within F' 
into channel k' G £q picking up a factor S'{E)j^/i. Then it propagates from T' to T" picking 
up a phase factor exp(z\/^afc'). To end up in channel k it then has to be transmitted within T" 
into another channel k" which is in £q and which we have identified with £q. Therefore it picks 
up a factor S"{E)kii and then a factor exp(z\/^Ofe//) when propagating back to V and finally 
comes the factor S'{E)k k" from propagation within T' before ending up in channel k. By the 
superposition principle summation has to be carried over all such k' and k". This contribution 
therefore corresponds to the term with m = in the expression for K2 in (32) when inserted 
into the second term in the first relation of (|^). The higher order contributions arise if the plane 
wave is reflected m + 1 > 1 times back and forth between V and T". Again by the superposition 
principle one finally has to sum over all m, giving the desired relation for S{E)ki. 

Example 4.1. Consider an arbitrary self-adjoint Laplacian A(A, B) with local boundary con- 
ditions on the graph depicted in Fig. ^ where the distance between the vertices is a. The 
composition rule ( 34) with 



V{a) 



easily gives 



Su 
S22 
S12 
S21 



11 "T '-'l2'-'ll'-'2lU ~ '-'22'-'ll6 ) 

^22 I '-'22'-'21'-'l2V ~ '-'22'-'ll^ ) 

Ql Qll l-\ Ql n// J2iay/E\ — l 
'-'l2'-'l2l-'- ~ '-'22'-'ll^ J ' 

qII qI (-t q/ q// J2ia\/^\ — l 
'-'21'-'2lU ~ '-'22'-'ll6 ) ) 



where the S-matrices are written in the form analogous to ( |30| ) 



Ql — ( ^11 ^12 1 qll 

\ Ql Ql ; ' 

V '-'21 '-'22 / 



Qll 

'-'21 



all 
'-'12 
an 
'-'22 



leaving out the E— dependence. These relations are equivalent to the Aktosun factorization 
formula applied to the Laplacian on a line with boundary conditions posed at x = and x = a. 



Example 4.2. Consider the Laplacian on the graph from Example 3.2 (Fig. 0j. The relation 



(34) holds with 



Via) 





V 





i\fEa n 
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Figure 3: The graph from Example 4.1. 
and S' = S" calculated in the Example 2.4. It is easy to check that 

Therefore from the formula f|7^ again follows. Note that Ki = K2 is singular at Tia,b> 
but these singularities disappear in the on-shell S-matrix. 



As already remarked multiple application of (34) to an arbitrary graph allows one by com- 
plete induction on the number of vertices to calculate its S-matrix from the S -matrices corre- 
sponding to single-vertex graphs. If these single vertex graphs contain no tadpoles, i.e. internal 



lines starting and ending at the vertex, then (|6|) and ( |34D give a complete explicit construction of 
the S-matrix in terms of the on-shell S-matrices discussed in Section 2. In case when a resulting 
single-vertex graph contains tadpoles we proceed as follows. Let the graph V have one vertex, n 
external lines and m tadpoles of lengths Oj. To calculate the S-matrix of T we insert an extra ver- 
tex on each of the internal lines (for definiteness, say, at x = aj/2). At these new vertices we im- 
pose trivial boundary conditions given by the choice a — l = d — l = & = c = 0, exp(2i;u) = 1 
of Example 3.2. With these new vertices we may now repeat our previous procedure. Thus 
in the end we arrive at graphs with one vertex only and no tadpoles. But the unitarity of the 
associated on-shell S-matrices was established in Section 2, so this property in the general case 
follows from Theorem 4.1. This is our third proof of unitarity. As an illustration we consider 
the following 




Figure 4: The graph F with n = 1 and m = 1: a closed loop (i.e. tadpole) of length a plus one 
external line. The open circle denotes the vertex added. 
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Example 4.3. Consider the Laplacian on the graph F depicted in Fig. ^with a circle of length 
a and the same boundary conditions at the 3 -vertex as in Example 2.4. Then 



S'tadpolo (-E) — e 



An easy calculation shows that 



St.dpoie{E) = S'i'^^E) *p=2 V{a)S{E)V{c 



where S{E) is the on-shell S-matrix of Example 2.4 and 

Via) 

In this case Ki equals 



( e'^^l'' 
e^^'^/^ 





1 



^ 3 



\fEa\-\ 



^v^a^-l ( 1 



2Jy/Ea 
3^ 

I iy/Ea 

3^ 




which is singular for E 



(27rn)^ 



The feature exhibited in this example (see also Example 3.2) may be generalized as follows. 
Let r be a graph with local boundary conditions at each vertex which are scale invariant in the 



sense of Corollary 2.3. Then by (34) the E dependence of the on-shell S-matrix St{E) for this 
graph only enters through the lengths a in the form exp{i\/^aj). In particular linia^o Sr{E) 
is independent of the energy. This observation might be helpful in deciding which boundary 
conditions might be physically realized in a given experimental context. 



5 Conclusions and Outlook 

In this article we have established unitaritity of the S-matrix for arbitrary finite quantum wires 
with a Hamiltonian given by an arbitrary self-adjoint extension of the Laplace operator. The ex- 
plicit determination of the on-shell S-matrix has been reduced to a finite matrix problem, which 
thus is accessible to computer calculations. A quantum wire with two open ends but arbitrary 
interior and arbitrary boundary conditions may be viewed as a theory with point interaction and 
an internal structure (see e.g. [ |85| , p^]). Ultimately relativistic, local quantum field theories pro- 
vide the appropriate set-up for considering point-like interactions and internal structures. Thus 
for example the <I>^-theory is the quantum field theoretic version of the (5-potential. The com- 
position rule established in Section 4 gives a generalization to arbitrary graphs of the Aktosun 



factorization formula [ 5C ] for potential scattering on the line. 



Our approach offers several generalizations. First it is possible to introduce a potential V 



on the entire wire thus replacing -A{A, B) by -A(A, B) + V (see e.g. [|8§], [|87|]). Although 
solutions for the on-shell S-matrix in closed form like above may not be obtainable in general, 
most structural properties should still hold, provided the potential is sufficiently strongly decay- 
ing at the infinities. Also one may construct lattice models in the following way. To each lattice 
site j G Z'^ in M'^ one may associate an S-matrix with 2k open ends and connect neighboring 



lattice sites accordingly (see e.g. [33, 3g, This would lead to the possibility of considering 



infinite quantum wires - obtained as 'thermodynamic' limits of finite quantum wires - and study 



their conductance properties in the spirit of e.g. the analysis in [33]. Thus if one takes the same 
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S-matrices at each site the resulting theory will be translation invariant. The generalized star 
product also offers the possibility of introducing a transfer matrix theory (see again [^]). For 



the case k = 1 this will generalize the Kronig-Penney model []34|,|89|]. Finally the S-matrices at 
the lattice sites could vary stochastically allowing the study of percolation effects. 

The present discussion of Laplace operators on graphs could also be used to study Brow- 
nian motion and the associated diffusion process. More precisely we expect the heat kernel 
exp(tA(A, B)){x, y) to be nonnegative if the boundary conditions {A, B) are real. It would be 
interesting to obtain a representation in the spirit of the Selberg and Gutzwiller trace formula. If 
the boundary conditions are not real, reflecting the presence of a magnetic field say, we expect 
a corresponding Ito formula to be valid and Aharonov-Bohm like effects to show up. 
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Appendix A 

In this Appendix we will describe A (A, B) defined on a graph with a single vertex from the 
viewpoint of von Neumann's extension theory (see e.g. [^]). According to von Neumann's 
theorem any self-adjoint extension /S.{A,B) of A'^ can be uniquely parametrized by a linear 
isometric isomorphism Wa,b '■ Ker(— A'^^ — i) ^ Ker(— A'^^ + i) according to the formula 

V{A{A, B)) = ^^p + ^+ + Wa,b^P+\ ^ G :P(A°), G Ker(-AO^ - i)} , 
-A{A, B){il, + ^+ + Wa,b^+) = -AV + - iWAM+- 

To express yVA,B in terms of the matrices A and B we compare A(y4, B) with another self- 
adjoint extension A(A',i?') of ISP. For definiteness we take A' = 1 and B' = 0, which 
corresponds to the Dirichlet boundary conditions (but another choice would also be possible). 
Let Uj G Ker(— A'^^ — i) and Vj G Ker(— A°^ + i), j = 1, . . . , n be given as 

(n,(x))fe=5,-fc2V4e75(-i+^)-, (z;,(x)), = 5,,2i/4e73(-i-^)-, k=l,...,n. (Al) 

One can easily verify that {u}"^^^ and {i'}j=i are orthonormal bases for Ker(— A°^ — i) and 

Ker(— A°^ + i), respectively. Denote by Wa,b the unitary matrix representation of Wa,b with 
respect to the bases and i.e. 

n 

y^A,BUj = ^ {WA,B)kj Vk- 
k=l 

It is easy to see that Wa'=i,b'=o = —I- Our strategy will be to apply Krein's formula (see e.g. 



||64|]) to obtain Wa,b in terms of the matrices A and B. To construct A{A, B) in the sense of 
von Neumann we proceed in two steps. In the first step let A be the maximal common part of 
A(^, B) and A{A' = I, = 0), i.e. A is the largest closed extension of A° with 

V(K) = V{A{A, B)) n V{A{A! = I,B' = 0)) = { V| V'(O) = 0, Bib'{0) = 0} . 
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Obviously, A is a symmetric operator with defect indices (r, r), r = Ranki? (0 < r < re). 
Also one has Ker(— A^^ =F ^) ^ Ker(— A*^^ =F i) such that we may decompose 

Ker(-A°"^Ti) = Ker(-AtTO©-^±- (A2) 

Performing, if necessary, a renumeration of the basis elements we may assume that {uYj^^ and 
{vYj^i are bases for Ker(— A"!" — i) and Ker(— A"!" + i) respectively. 

Von Neumann's theorem describes A as a symmetric extension of A*^. All such extensions 
are parametrized by linear partial isometrics W' : Ker(— A*^^ — i) — > Ker(— A"^ + i). For 
the case at hand it is easy to see that the isometry corresponding to A defines an isometric 
isomorphism W: Ai+ ^ Ai- and therefore 

P(A) = {v- + -0+ + VVV'+I V' G ^?(A°), i^+eM+Y 
-A(V' + tl'+ + yV^p+) = -A^V + itp+ - iyV'4'+- 

Denoting by W the matrix representation of W in the bases {n}"^^_|_^ and {v}'j^^_^_^ we get that 

W = -I. _ 

In the second step let Wa.b '■ Ker(— A^ — i) — > Ker(— A^ + i) be the linear isometric 
isomorphism parametrizing A{A,B) as a self-adjoint extension of A. Denote by Wa,b its 
unitary matrix representation with respect to the bases {uYj^^ and {vYj=i- Again it is easy to 

see that Wa'=i,b'=o = — I (r x r matrix). Also from the discussion above it follows that 

Wa,b = Wa,b © (-1) (A3) 



with respect to the decomposition ( |A2| ). 

Krein's formula (see e.g. ||6^) now states that the difference of the resolvents of —A{A, B) 
and —/S.{A' = 1,B' = 0) at the point z = i is given by 

r 

RA,B{i) - Ra'=i,b'=o{'i) = Pjkivk,-)vj, (A4) 

j,k=i 



with some r x r matrix P of maximal rank. By Corollary B.3 of [ |90| ] it follows that 

P='-{I + W^i,). (A5) 



Obviously, relation ( |A4| ) can be rewritten in the form 

n 

RA,B{i) - RA'=I,B'=Q{i) = Pjk{vk,-)vj, (A6) 

j,k=l 

with P being n x n matrix of rank r such that P = P © with respect to the decomposition 
( |A2| ). Thus from (Q and (Q it follows that 

p='-{i + wxi,). (^"7) 

The resolvent RA'=i^B'=o{i) of A{A' = l,B' = 0) can be given explicitly. Its integral 
kernel (Green's function) for x < y has the form 

RA<=i,B'=Q{x,y;i)=\ T^e^^ '\ 
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such that 



RA'=i,B'=o{0,y;i) = 0, 



ORa'. 



'B'=0 



dx 



(0,y;i) 



Since the boundary conditions take the form 



araM^^ y; + = 



for all y > in terms of the resolvents, from (A6) we obtain 



1 



A + ^(-1 + i)B]P = — -B 



V2 



1 



^/2 



Comparing this with (A7) we obtain 



W 



A,B 



V2' 



-1 



A + ^{-l + i)B] {A-^{l + i)B 



V2' 



Note that {A, B) given as 



1 - 1 
A ^B, B = -^B 

V2 \/2 



also defines a maximal isotropic subspace and VF^^ = g{E = 1) (the S-matrix for the sin- 
gle vertex theory) holds, such that in particular Wa^b is unitary and satisfies Wa,b = WcA,CB 
for any invertible C, as it should be. As an example, for A = and B = 1 (Neumann boundary 
conditions) this gives Wa=o b=i = again as it should be. 



Appendix B 

Here we give an alternative 'analytic' proof of Theorem 3.3, i.e. of the unitarity of the on-shell 
S-matrix. We resort to an argument which is a modification of well known arguments used 
to prove orthogonality relations for improper eigenfunctions of Schrodinger Hamiltonians (see 
also the remark at the end of this appendix) and which is instructive in its own right. For this 
purpose we introduce Hilbert spaces H r indexed by i? > and given as 

where the spaces Tii are as before and where He,R = i^([0, R]) for all e&£. Then E) E 
Hr for all e G <5 and all < i? < cxd. Let { , )r denote the canonical scalar product on Hr. 
Also we have Aijj^{-, E) = -E^^{-, E) such that for all E,E' Y.a,b and all k,l e £ 

{A^\;E),^\., E'))r - ii^'^i; E), Ai;\.,E'))R 

= {E' -E){^1;\.,E),^\.,E'))r. (B1) 

Now all terms in this equation are smooth in £" G M \ Tia,b and R. We may therefore take 



the limit E' E. This gives for the right hand side of ( |Bl| ) 



lim - e){'4j\;E),^\-,e'))r = 0. 

E — >Ej 
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On the other hand we may evaluate the left hand side of by performing a partial integration. 
Since both E) and il/{-,E') satisfy the boundary conditions {A, B), the only contributions 
arise at the n 'new' boundaries x = R, i.e. 

(AV^^(-, E),^\; E'))r - E), A^'(-, E'))r 

ee£' 

- ^)e'i^+^)^ Ski{E'). 
The limit of the right hand side as E' ^ E equals 

Vee£' / 

which by the preceding arguments is zero, concluding this second proof of unitarity. 

We could have chosen E' = E from the very beginning. However, the present discus- 
sion may be applied to show general orthogonality properties for the improper eigenfunctions 
ip^{-,E). This is achieved by first dividing ( pT| ) by E' — E and then taking the limit i? — > oo in 
the sense of distributions in E' and E. 



Appendix C 

This appendix is devoted to a proof of Theorem 4.1. We start with the following observation. 
For U' written in block form as in ( |30| ) we define a map U' U'^ with r = T{n' — p,p) given 
as 



U' U' 

^22 ^21 

U[2 U[, 



which amounts to interchanging the first n' —p indices with the last p indices while keeping the 
order of the indices otherwise fixed. U"^ and U"^ are defined analogously with r = r(p, n" —p) 
and r = r(n' — p, n" — p) respectively. From the definition of *v it follows immediately that 
the following 'transposition law' holds 

[/- = U"^ *v-i U'^ 

whenever U = U' *v U". To prove Theorem 4.1 we have to show that the following four 
relations hold 

ul,Un + ulU2i = I, 

UlUu + UlU22 = 0, 

UI2U12 + UI2U22 = I, 
ul^Un + ul^U2i = 0. 
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By our previous observation it suffices to prove only the first two relations. To prove the first 
one, we insert the definition of U and obtain 

ulUn + ulU2i = U[\U[^ + U[\U[2K2U'^^VU^^ 

+u!,{v-'u';Ik2U[\u[, 

+u!,\ v-'u';i kI U[\ U[2 K2 t/fi V [7^1 

5 

= IZ"- 

i=l 

We now use the unitarity relations for U' and U" . This gives 

5 

i=2 

+V-^ U'll kI [7^2 U[2 K2 U^i V + k\ U'^ U'^^ Ki}U^^ 
= U'l{-U!22 K2 V - U'll kI + V-^ U'll K2 C/fi V 
-V~^ U'll kI u!^2 U^2 K2 U^i V + k\ Ki - k\ U';l U'{^ Ki}U^^. 

To establish Ylf=i «i = I it therefore suffices to show that the expression in braces equals I. 
Using one of the relations in and its adjoint we have 

kIKi = l + U22K2U'{-^V + V-^U'(IkIu'^^ 
+v~^ u'll kI U'^^ U22 K2 U'{^ V 

and 

Klu';lu^,K^ = y-i(i + c/('ti^t[/^ty-i) 

U';lu'{^{l + VU'22K2Ul)V. 

Hence it suffices to show that 

-(I + u'll kI f/22 v-^)u'll u'l^{i + vu'22 K2 u'l{) + u"l kI K2 u'l^ = 0. 

To show this it suffices in turn to prove that 

kIK2 = 1 + Klu'^^V-^Ull + U'{^VU22K2 

+kI u'^^ v^^ u'll U'l^ V U22 K2. 

But this relation follows by inserting the relation 

K2 = y-^ + y-^ Uu y U22 K2 

and its adjoint into the left hand side. To sum up we have proved the first of the unitarity relations 
for U. To prove the second relation we again insert the definition of U, use the unitarity relations 
for U' and U" twice and obtain 

ul^Ui2 + ul^U22 = U!,\{-U^2K2 + y-^U'lKlK2 

-y-^u';iKlu'lu'22K2 
-k\u'^ + k\kiU'22 y-^ 
-Klu';iu^^KiU'22y-^]u^2- 
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Hence it suffices to show that the expression in braces {= a^) vanishes. But 

-kI U';l - kI U'^ U';^ Ki = -kI U'll V K2 

= -v-\i + u"l kI u^^ v-^)u"I V K2 

= -V-^ U'^l VK2- V-^ U'll {kI - V)K2 

= -v-^u"IkIk2 

implies 

ae = -U^2 K2 - V-^ C/(;^ kI U'^^ U'22 K2 + k\ U'22 V'K 
The chain of equaUties 

k\ k, U'22 v-^ = + u';i kI u'l){i + u'22 K2 u^i v)u'22 v-^ 

= U'22V-' + U'22K2U'{^VU'22V-^ 

+v-^ u'-ii kI [/22 u'22 K2 w^i V u'22 

+v-'u';IkIu^\u^^ 

= U'22V-^ + U^2{K2'V-^) 

u'^IkI u'l u^2{K2 - v-^) 
+v-'u';IkIu!^^u^2V-' 

= U'22K2 + V-^U';IkIu'^^U'22K2 

finally leads to ae = as desired concluding the proof of Theorem 4.1. 
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